Long term ordering kinetics of the two dimensional q-state Potts model 
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We studied the non-equilibrium dynamics of the q-state Potts model in the square lattice, after 
a quench to sub-critical temperatures. By means of a continuous time Monte Carlo algorithm 
(non-conserved order parameter dynamics) we analyzed the long term behavior of the energy and 
relaxation time for a wide range of quench temperatures and system sizes. For ^ > 4 we found the 
existence of different dynamical regimes, according to quench temperature range. At low (but finite) 
temperatures and very long times the Lifshitz-Allen-Cahn domain growth behavior is interrupted 
with finite probability when the system stuck in highly symmetric non-equilibrium metastable states, 
which induce activation in the domain growth, in agreement with early predictions of Lifshitz [JETP 
42, 1354 (1962)]. Moreover, if the temperature is very low, the system always gets stuck at short 
times in a highly disordered metastable states with finite life time, which have been recently identified 
as glassy states. The finite size scaling properties of the different relaxation times involved, as well 
as their temperature dependency are analyzed in detail. 
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I. INTRODUCTION 

The problem of domain growth kinetics in systems with 
many degenerate ground states had attracted a lot of at- 
tention in the past^s^J^i^i^i^. In early works Lifshitz^ and 
then Safran^ posed the discussion about the effect of ac- 
tivated processes in the domain growth. They suggested 
that d-dimensional, g'-state degenerate models could be- 
come trapped in local metastable states for g > d + 1, 
which would then greatly slow down the relaxational ki- 
netics. Their argument in two dimensions is that a hon- 
eycomb structure of hexagonal domains is stable under 
small distortions of the interfaces, because such distor- 
tions do not increase the interfacial free energy. Large 
distortions due to fluctuations are then needed to move 
the interfaces and so domain growth becomes activated. 
The prototype of such a system is the q-state Potts 
model^, whose Hamiltonian is given by 



H = -J ^5{si,Sj) J>0 



(1) 



where si = 1, 2, . . . , g', 5{si^ Sj) is the Kronecker delta and 
the sum runs over all the pairs of nearest-neighbor sites. 
Vinals and Grant^ performed Monte Carlo (MC) simula- 
tions with Glauber dynamics for g = 8 in square lattices 
of size up to = 128^ sites and time scales of the or- 
der 10^ Monte Carlo Steps (MCS) (1 MCS is defined as 
as a complete cycle of N spin update trials, according 
to a Metropolis MC algorithm). They found evidences 
of metastable configurations composed by squares of dif- 
ferent colors (i.e., different values of q) that freeze the 
dynamics a T = 0. Moreover, they argue that those 
configurations only effectively dominate at T = and 
growth at finite temperatures is not limited by activated 
processes^. Hence, for long enough times the average 



linear domain size should follow the Lifshitz-Allen-Cahn 
(LAC) law l{t) ~ t-^/^. Grest and coworkers^ reported 
results using MC simulations with N = 1000^ sites and 
values of g up to g = 64 at very low temperatures, which 
are consistent with LAC behavior for any value of q, at 
least for time scales up to 10^ MCS. Analytical results on 
a coarse grained model also confirm those results^. How- 
ever, the time scales considered in those works are very 
short to exclude the existence of activated process of the 
type predicted by Lifshit^ and Safran^. Moreover, re- 
cent investigations^ii^ about the pinning configurations 
found by Vinals and Grant^ showed that the system gets 
stuck in those states at times scales longer than 10^ MCS; 
pinning also happens at low but finite temperature^, 
although in that case the metastable states present a 
finite life time that increases with— q. Those works 
also showed that the nature of those highly disordered 
metastable states is more related to a glassy ono^i^^^ 
than to the type of configurations predicted by Lifshit^ 
and Safran^. We will refer hereafter to those disordered 
metastable states as the glassy ones. 

In this work we concentrated mainly in the q = 9 
case in square lattices with N = L x L sites and peri- 
odic boundary conditions. Some complementary results 
are also presented for other values of q. The imple- 
mentation of a continuous time MC algorithm or n-fold 
technique^^'^^ allowed us to analyze a large statistics on 
system sizes up to L = 500 and time scales running from 
10^ to 10^^ MCS. Some details of the continuous time 
MC implementation are presented in Section [III 

We analyzed the relaxation of the system after a 
quench from infinite temperature down to subcritical 
temperatures T, for a wide range of values of T and L. 
We found the existence of different dynamical regimes, 
according to the quench temperature range. First of all 
there is some characteristic temperature T* < Tc, such 
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that for T* < T < Tc simple coarsening dominates the 
relaxation (that is, the domain growth follows the LAC 
law even for very long times), except very close to Tc, 
where evidences of nucleation relaxation mechanisms ap- 
pear; however, the last case will not be analyzed in this 
work and the corresponding results will be presented in 
a forthcoming publication. For T < T* we found that, 
at long time scales (i.e., much longer than those con- 
sidered in previous works), the LAC relaxation is inter- 
rupted when the system gets trapped in highly symmetric 
metastable states with finite probability (i.e., for a large 
fraction of realizations of the thermal noise that do not 
decrease as the system size increases) . We found two dif- 
ferent types of those configurations: striped states and 
honeycomb like structures; the latter are configurations 
composed by macroscopic six-sided irregular polygonal 
domains of different colors. Striped states are composed 
by two macroscopic ferromagnetic domains with straight 
walls parallel to the coordinate axis and has been pre- 
viously observed in the Ising mode]~^^^^^ {q = 2) 
at T = 0. For (7 > 3 and T = the probability to 
reach a striped state becomes zero in the thermodynamic 
limit^^; we found that low but finite temperatures make 
that probability to become finite. The presence of hon- 
eycomb like structures in the q-state Potts model, as far 
as we know, has not been previously reported and is in 
agreement with the Lifshitz's prediction^. When the sys- 
tem reaches either striped or honeycomb like states the 
dynamics becomes activated. 

Finally, we found a temperature <C Tc, such that 
for T < Tg the system always gets stuck at inter- 
mediate times in a glassy state of the type reported 
previously^i^ii^iii for other values of g > 4. We found 
that for q = 9 those states present a finite life time (i.e., 
independent of the system size) with a well defined free 
energy barrier associated to it. 

We also analyzed the scaling properties of the char- 
acteristic times associated with the different relaxation 
processes, as well as the probability of reaching a striped 
or honeycomb state for large values of q. 

We verified that the whole relaxation scenario is quali- 
tatively observed when open(instead of periodic) bound- 
ary conditions are used. All the numerical results are 
presented in Section [1111 and some discussion is presented 
in Section HVl 

Besides its theoretical interest, the large-g Potts model 
(or variations of it) is used for simulating the dynamics of 
a large variety of systems, such as soap froth^^'^^, grain 
growth^^'^^ and biological cells24. The present results 
help to establish the conditions under which equilibrium 
can be actually reached, as well as the different possible 
low temperature relaxation scenarios. 



described by Hamiltonian ([T]), after a quench from infi- 
nite temperature (i.e., a completely random initial con- 
figuration) down to sub-critical temperatures. The Potts 
model undergoes a second order phase transition for 
g = 2, 3, 4 and a first order one for q > 4:, where the crit- 
ical temperature in the square lattice is known exactly^ 
for any value of q and is given by Tc = 1/ln (l + y^) 
(hereafter we will use natural units ks = J = l)- Most 
of our analysis were concentrated in the q = 9 case, for 
which Tc = 0.7213 .. ., and some complementary calcula- 
tions were performed for g = 2, 3, 4, 5, 15 and 30. 

We were mainly interested in the late stages of the dy- 
namics, where large domains are formed. In that case, 
the computational cost of a single spin fiip dynamics (for 
instance, heat bath algorithm) becomes very high, be- 
cause the fiipping probability of spins inside the domains 
(which are the majority) becomes very small. An efficient 
way to achieve longer time scales for reasonable large sys- 
tem sizes is the usage of continuous time MC methods 
or n-fold techniques'^. In these type of algorithms a 
fiip occurs a each step and the time such event would 
have elapsed in a single spin fiip algorithm is calculated 
from the associated fiipping probability. Let us briefiy 
summarize the implementation of the algorithm for the 
nearest-neighbors g-state Potts model. For a given spin 
configuration we will call "potential spins" to the q — I 
possible states for each site in the lattice, different from 
the present ones. All the potential spin of the system 
are then classified into lists, where members of a given 
list would produce the same change in the energy of the 
system, if chosen to replace the old spin state in the cor- 
responding site. For a single fiip there exist only 9 pos- 
sibilities (^ = -4,-3,-2,-1,0,1,2,3,4), so we have 
9 classes for any value of q > 3. In the Metropolis algo- 
rithm the probability of a potential spin belonging to the 
class / (/ = 1, . . . , 9) to be effectively fiipped is: 



Pi 



ni 



N{q-1) 



mm 



1, exp 



kBTj 



(2) 



where ni is the number of spins in class /. The total 
probability of any fiipping event occurring in a given step 
is 



(3) 



In the present algorithm, at each step a class is sorted 
with probability ([2]) and a potential spin is sorted with 
equal probability among all the members in the class. 
After updating the corresponding spin and the lists, the 
time step is incremented by an amount 



II. METHODS At = — ^ U) 

NQ ^ ^ 

We analyzed the time evolution under a type A dy- where r is a random number uniformly distributed be- 
namics (non-conserved order parameter) of the system tween zero and one and the time step At is measured in 
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FIG. 1: (Color on-line) Relaxation function 
^ = 9 and temperatures ranging from T - 
T = 0.72 (top). 



6£;(t)forL = 300, 
: 0.5 (bottom) to 



MCS. The details of the algorithm can be seen in Ref^^^. 
The implementation of this algorithm allowed us to per- 
form simulations for time scales running from 10^ MCS 
(for sizes up to L = 500) to 10^^ MCS (for sizes up to 
L = 100). In order to check the algorithm we also re- 
peated several of the simulations using a single spin-flip 
algorithm (heat bath) for L = 200 and time scales up to 
10^ MCS. The results were identical. 

Our analysis of the dynamics was mainly focused on 
the behavior of two quantities: the average energy per 
spin e{t) = {H{t)) /N as a function of t (the average was 
taken over different initial conflgurations and different re- 
alizations of the thermal noise) and the equilibration time 
r; the last quantity was deflned as the time at which the 
instantaneous energy falls below an equilibration thresh- 
old. Such threshold was set as the equilibrium energy at 
the corresponding temperature plus one standard devia- 
tion, where those quantities were flrst calculated by run- 
ning a set of simulations starting from the ordered state 
and letting the system to equilibrate. We calculated the 
probability distribution (normalized histogram) P{r) for 
different values of T and L. 



III. RESULTS 

In order to compare the behavior of the average en- 
ergy per spin e{t) = (H) /N for different temperatures, 
we flrst introduce the relaxation function (or normalized 
excess of energy) 



e{t) — e{oo) 
e(0) -e(oo) 



(5) 



where e(oo) is the equilibrium energy. In FiglDwe show 
the typical behavior of (pEit) for L = 300 and tempera- 
tures between 0.72 and 0.5 (T^ = 0.7213 . . .). For temper- 
atures close enough to Tc (0.715 < T < Tc) we see that 
the system clearly stuck in a high energy metaestable 
state. Close examination of different quantities in the 
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FIG. 2: (Color on-line) Equilibration time probability distri- 
bution P{r) for L = 100 and ^ = 9. (a) Temperatures ranging 
from T = 0.4 (top) to T = 0.719 (bottom); (b) temperatures 
ranging from T = 0.3 (left) to T = 0.09 (right). 



metastable state show that this corresponds to a disor- 
dered (i.e., paramagnetic) one and hence it is directly 
related to the flrst order nature of the transition. More- 
over, we found evidences that in this regime relaxation 
is dominated by nucleation mechanisms, but the details 
of that analysis will be presented in a forthcoming pub- 
lication. For temperatures T < 0.715 we see that the 
metastable plateau disappears and the relaxation func- 
tion decays (after a short transient) for all temperatures 
as ^£;(t) ~ Since the excess of energy respect to 

the equilibrium state in a domain growth process is given 
by the average energy of the domain walls, a simple cal- 
culation shows that 0£;(t) ~ ^/l{t), l{t) being the average 
linear domain size. Hence, the behavior of FigHJis consis- 
tent with the LAC law. As we will show later, the flnite 
size scaling properties of the average typical equilibration 
time in this temperature range are also consistent with 
the LAC law. 



A. Relaxation at intermediate temperatures and 
blocked states: characterization and scaling 

In Figl2] we see the typical behavior of P{r) for an 
intermediate size {L = 100) and different temperature 
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FIG. 3: (Color on-line) Characteristic relaxation times n and 
T2 vs. 1/T for L = 100 and ^ = 9; the same qualitative 
behaviors are observed for other system sizes. The continuous 
lines are a guide to the eye. 
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FIG. 4: (Color on-line) Energy per spin as a function of 
time and typical spin configurations in one realization of the 
stochastic noise, when the system gets stuck in a striped con- 
figuration (L = 200, ^ = 9 and T = 0.2). The different colors 
codify different spin values = 1, . . . , 9 



ranges. The different dynamical regimes can be already 
appreciated in this figure. Close enough to Tc (T = 0.719 
in Figl2^) P(r) exhibit a well defined peak centered at 
a characteristic time Tnuci ^ 10^ MCS, which is associ- 
ated to a nucleation based relaxation mechanism already 
mentioned. As the temperature decreases below some 
temperature 0.715 < Tn < Tc, this peak is suddenly re- 
placed by another one centered at a characteristic value 
Ti , which is about one order of magnitude smaller of Tnuci 
and remains almost independent of the temperature in 
the range 0.3 < T < 0.715; in the temperature range 
< T < 0.2 (see FiglSt), ri exhibit a strong temperature 
dependency. For temperatures smaller than (but close 
to) Tn (see Figl2^), P(r) develops a long right tail; for 
temperatures T < T* ^ 0.6 the tail becomes a distinct 
peak centered at a new characteristic time r2, which in- 
creases exponentially as the temperature decreases. This 
behavior indicates the existence of two distinct phenom- 
ena affecting the relaxation at different time scales, where 
T* acts as a reference temperature signaling the time 
scales separation crossover point. The temperature be- 
havior of Ti and r2 is summarized in the Arrhenius plot 
of Figl. 

We will show that n is associated with simple coars- 
ening processes that follow LAC law for all times, while 
r2 is associated to processes in which the system gets 
stuck in striped met ast able states, composed by two fer- 
romagnetic states whose walls are parallel to coordinate 
axis, as shown in the example of FiglH Those type of 
metastable states have been already observed in the two 
dimensional Ising model {q = 2) at zero temperature, 
where they become fiozen^^^i^^. At finite tempera- 
ture, striped states perform a random parallel movement 
in the direction perpendicular to the walls. Hence, in a 
finite system those states relax to equilibrium when both 
walls collapse. Spirin, Kaprivsky and Redner— showed 
that the basic mechanism for the parallel movement of a 
straight domain wall is the creation of a "dent", that is. 



the flip of one of the spins adjacent to the wall. Since 
after flipping the spin its neighbors can flip without en- 
ergy cost, the energy barrier for the creation of a dent 
is 2 (in units of the coupling constant J) for the q = 2 
Potts model (or 4 for the Ising model). For q > 2 the 
energy cost of any other movement (including a flipping 
to a third color different from those of the domains) is 
larger. Hence, once the striped state is reached, the time 
needed to relax should be basically independent of q and 
this is consistent with the Arrhenius behavior r2 ~ e^/^ 
observed in Figl3l 

From Figl3]we can notice also that, for a wide range of 
temperatures T < Tn (approximately down to T ^ 0.2) 
Ti remains almost independent of T, consistently with 
a simple coarsening behavior; at lower temperatures we 
see a crossover into an activated behavior, that will be 
analyzed later. 

A deeper understanding of the mechanisms involved in 
the relaxation can be obtained from the finite size scaling 
of the different quantities involved. In Figl5]we show the 
typical behavior of P{r) for different system sizes at a 
fixed temperature T = 0.2. The first thing we note is 
that the two-peak structure remains in the large L limit. 
Moreover, the ratio between the areas below both peaks 
becomes constant in such limit. The same property is 
observed for temperatures up to T*. We will analyze 
this in more detail at the end of this section. Let us now 
consider the finite size scaling of the relaxation times. 

From FiglG] we see that n ~ for a wide range of 
temperatures, both above and below T*. This is also con- 
sistent with a simple coarsening process, in which equili- 
bration will be attained once 

Let us now analyze the finite size scaling of r2. Spirin, 
Kaprivsky and Redner^^ suggested that, at low enough 
temperatures, the movement of a flat interface will be 
dominated by processes involving a single dent creation; 
once the dent is created it performs a random walk, un- 
til either the dent disappears or it covers the whole line, 
where the typical time needed for the last event scales 
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FIG. 5: (Color on-line) Equilibration time probability distri- 
bution P(t) for T = 0.2, q = 9 and system sizes ranging from 
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FIG. 7: (Color on-line) Characteristic relaxation time T2 vs 
L for T = 0.2 and different values of q. The dashed lines 





FIG. 6: (Color on-line) Characteristic relaxation time n vs L 
for Q = 9 and different temperatures; dotted lines are a guide 
to the eye. 



FIG. 8: (Color on-line) Finite size scaling exponent cu of the 
characteristic relaxation time r2 as a function of T for different 
values of q; dotted lines are a guide to the eye. 



as^^ L. This mechanism leads to a random walk move- 
ment of both interfaces, so there must be typically 
such hopping events for the interfaces to meet and there- 
fore the relaxation time should scale as^^ L^. However, 
this argument only works for small system sizes. Once a 
dent is created, the probability of creation of new dents 
along the interface, before the dent covers the line, in- 
creases with the system size; therefore, the typical time 
for a one-site hoping event of an entire interface should 
increase slower than linearly with L and r2 slower than 
L^. This can be appreciated in the clear crossover from 
l3 to r2 ~ with a; < 3 around L = 15, ob- 
served in FiglJl both for g = 2 and q = 9 (the same 
effect is observed for any temperature T < 0.2). Striped 
states appear for temperatures up to T*, but the the walls 
movements are no longer dominated by one-site hopping 
events for T > 0.2; instead of that, a direct inspection of 
the spin configuration during relaxation shows that, at 
temperatures close to T* the movements of the domain 
walls resemble (for large system sizes) that of elastic lines 
subjected to a random noise. Hence, the temperature de- 
pendence of r2 departs from the e^/^ behavior, as can be 
seen from FiglH However, the finite size scaling r2 ~ 
still holds for temperatures up to T*, where the expo- 



nent uj displays a marked increase with the temperature, 
reaching values slightly larger than 3 as T approaches 
T* (see FiglH]). Those values of the exponent can be un- 
derstood through the following argument. Suppose that 
each line behaves as a chain of L unit masses joined by 
springs, constrained to move along the direction perpen- 
dicular to the wall and subjected to independent white 
noise. By solving the corresponding Langevin equations 
in the overdamped limit, a simple calculation shows that 
the distance between the centers of mass of both chains 
performs a Brownian motion with an effective diffusion 
coefficient that scales 3.8 D ^ . Since the distance be- 
tween walls is of the order of L, this implies that the typ- 
ical time needed for the walls to encounter should scale 
approximately as L^. For q = 2 Lipowski^^ has shown 
that this scaling holds even for relatively large values of 
the temperature T/Tc(2) ^ 0.8. The results of FigElsug- 
gest that the scaling properties of r2 are independent of 

showing that large degeneracies in the ground state 
have no influence in this relaxation process. 

Let us return to the equilibration time probability dis- 
tribution P(r). Another salient feature of this distri- 
bution for temperatures T < T* is that the right peak 
broadens for large system sizes. To show this we re- 
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FIG. 9: (Color on-line) Equilibration time probability distri- 
bution P{r) starting from different initial configurations (ar- 
bitrary normalization) for L = 300, T = 0.2. The inset images 
show typical blocked spin configurations at the corresponding 
times. 



draw P(r) for L = 300 and T = 0.2 in Figl (full line). 
A careful inspection of individual processes shows that, 
such broad peak is actually associated with two different 
types of met ast able configurations: the striped ones al- 
ready described and honeycomb like structures; the latter 
are composed macroscopic six-sided irregular polygonal 
domains of different colors (see inset of FiglS]), where the 
angles between domain walls at the three-fold edges fluc- 
tuate around 120^. Those states are in agreement with 
Lifshitz prediction^ for q > 3 and we shall call them Lif- 
shitz states. By a calculation of the equilibration time 
starting directly from the Lifshitz and striped states, we 
verified that the broad peak of P{r) is actually a superpo- 
sition of two peaks, each one with its own distinctive max- 
imum at characteristic times r2, for the striped states, 
and Ta for the Lifshitz ones (see FiglS]). Lifshitz states 
are only detectable for system sizes L > 100. Actually, 
isolated three-fold vertex between flat domain walls of the 
type predicted by Lifshitz^ also appear for smaller sys- 
tem sizes, but complete honeycomb-like structures can be 
stabilized during detectable time scales (i.e., larger than 
the characteristic coarsening time scales) only for large 
enough system sizes. To determine the scaling properties 
of Ts, we calculated the escape time probability distribu- 
tion starting from the closest configuration to a Lifshitz 
state, that is, from an almost perfect four-colored hon- 
eycomb configuration (commensurability with the system 
size does not always allow a perfect honeycomb structure) 
for different values of L and T; we show an example in 
Figini We verified that the system quickly relaxes from 
that configuration into a Lifshitz state, from which it can 
either relax directly to the equilibrium state or pass first 
to a striped state, giving rise to a second peak in the 
corresponding probability distribution (see FiglS]). For 
completeness, we also calculated the escape time prob- 
ability distribution starting from a perfect two-domains 
striped state; the result is also shown in FiglSl 

The temperature dependency of ts is shown in FigHOl 
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FIG. 10: Characteristic relaxation time rs vs temperature for 
^ = 9 and L — 100; the continuous line is a guide to the eye. 
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FIG. 11: (Color on-line) Basic relaxation mechanism of a Lif- 
shitz state at low temperatures. The circle in (b) marks the 
creation of a dent in the vertex of (a), by flipping a spin from 
2^3. Circles in (c) and (d) exemplify spins that can flip 
without energy cost along the three converging walls. The 
whole process may lead to the hopping of the whole struc- 
ture (vertex plus walls) and ultimately to the collapse of two 
vertex. 



where we see that it displays a clear Arrhenius behav- 
ior with an activation barrier of height one, which is the 
minimum possible energy barrier associated with a single 
spin flip. This can be understood if we analyze the basic 
mechanisms behind the relaxation from Lifshitz states. 
We observed that Lifshitz states relax when two vertex 
of an hexagon edge collapse. A vertex movement, with 
the consequent displacement of the converging walls, oc- 
curs through a series of random hopping events. In FiglTTl 
we show an example of the hopping of a vertex one site 
to the right; one-site hopping events in the other direc- 
tions follow a similar process with the same energy cost. 
The movement of a vertex starts with the creation of a 
dent, by flipping one of the spins located at the neighbor 
sites of the vertex, as depicted in FiglTTb. This move- 



ment has an energy cost of one unit. Once the dent is 
created, the neighbor spins at the three converging wahs 
are free to flip without energy cost (see FigsITTh andfTTH). 
generating a diffusive motion of the dent along the three 
hues, and may lead to the displacement of the whole 
lines. This hopping movement of the vertex ultimately 
leads to the collapse of two of them and the consequent 
disappearance of the Lifshitz state. The whole mecha- 
nism is completely similar to that described by Spirin, 
Kaprivsky and Redner^^ in the case of a flat wall be- 
tween two striped domains, except that the creation of 
a dent adjacent to a vertex has an energy cost of just 
one energy unit (instead of 2, as in the case of a dent in 
a flat interface), which explains the behavior of FiglTOl 
Since hexagonal domains in a Lifshitz state are macro- 
scopic, the same flnite size scaling arguments used by 
Spirin, Kaprivsky and Redner^^ for the relaxation time 
apply in this case. Hence, one expect rs ~ e^/^. For 
the system sizes available, we verifled this scaling at low 
temperatures with an exponent c<; ~ 3, but we would ex- 
pect this value to be reduced for larger system sizes, as 
in the case of striped states (r2). 



B. Probability of blocked states 

We analized the probability Pb{q) of getting stuck in a 
blocked state. We deflned blocked states as those char- 
acterized by flat walls between domains. For q > 3 this 
includes Lifshitz and striped states. For q = 2 the system 
can also get trapped in another type of blocked states, 
characterized by diagonal stripes, whose interfaces fluc- 
tuate without energy cost^''; we shall call them diagonal 
states. For q > 3 we did not observed diagonal states at 
flnite temperatures. Although their presence for q > 3 
with low probability cannot be excluded, probably they 
are replaced at flnite temperature by the Lifshitz states. 

From the previous calculations of P{r) we could esti- 
mate Pb{q) by deflning for (for every value of T and L) 
a threshold value rt{T^L), such that a single realization 
with r > is attributed to the presence of a blocked 
state; rt can be estimated as the flrst minimum of P{r) 
located above n (see for instance FiglSJ. This procedure 
reduces the calculation of Pb{q) to a binomial experi- 
ment, hence, a simple calculation shows that a sample 
size of 2000 runs is enough to guarantee a statistical error 
smaller than 1% in all cases, thus saving a lot of CPU 
time. 

In FigIT2] we show the results for g = 9. The main 
source of error in this calculation is the choice of rt , which 
is not always evident, due to large fluctuations in the 
histograms for small sizes and very low temperatures; 
the error bars in FigIT2l were estimated by varying rt. 
From FigIT2k we see that, at T > 0, P5(9) saturates in a 
flnite value for L > 100, indicating a flnite probability in 
the limit L ^ oo. In FigIT2b we show the temperature 
dependency of the saturation value. We see that A (9) 
goes to zero as T ^ 0, consistently with the results of 
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FIG. 12: (Color online) Probability of getting stuck in a 
blocked state Pb{q) for q = 9: (a) as a function of L for dif- 
ferent temperatures; (b) as a function of T for different sizes. 
The dashed line corresponds to a linear fitting of the points 
for L = 200, giving an extrapolated value of 0.008 di 0.01 at 
T = 0. 



Spirin, Kaprivsky and Redner^^. 

Next we analyzed the probability Pb{q) as a function 
of q. The results are shown in FiglT3]for T = 0.15 and 
values of q ranging from g = 2 up to = 30 . For q = 2 
and T = the probability of reaching a striped state 
•^j26 2/3^ while the probability of reaching a diagonal 
state i^ ^ 0.04. At T = 0.15 we found the values 
^ 0.345 and ^ 0.045 respectively, giving rise to the value 
Pb{2) ^ 0.39. The diflFerences with the T = values 
are consistent with the enhancement of the probability 
at flnite temperature, already observed for q = 9. 

For q>3 the probability Pb{q) falls down to a temper- 
ature dependent flnite value, that is almost independent 
of q and smaller than half of Pb{2). 

It is worth noting that Spirin, Kaprivsky and Redner— 
reported another type of blocked states for = 3, charac- 
terized by both straight walls and diagonal walls, the lat- 
ter fluctuating without energy cost; they call these states 
"blinkers". We did not observe blinkers at flnite tempera- 
ture, at least for periodic boundary conditions. Although 
their existence with low probability cannot be excluded, 
probably they decay into Lifshitz states in time scales 
smaller than the characteristic Lifshitz relaxation times. 



q 

FIG. 13: Probability of getting stuck in a blocked state Pb{q) 
for L 200 and T = 0.15. 



C. Low temperatures relaxation: glassy states 

Let us now analyze the coarsening at very low tem- 
peratures. The increase of n for temperatures T < 0.2 
observed in Fig 13] indicates that the normal coarsening is 
affected by some kind of activated process. The increase 
in this relaxation time is associated with the plateau dis- 
played by the relaxation function in FigHH We found 
that this plateau appears below some characteristic tem- 
perature 0.1 < < 0.2 for q = 9. This plateau corre- 
sponds to a disordered metastable state characterized by 
almost square-shaped domains with a wide distribution 
of sizes (see inset in Fig IT5|) . That type of metastable 
state was previously reported for q = 7 and it was iden- 
tified as a glassy one^'^'^^. These states are only present 
for— q > 4. We verified that for g = 9 the normal coars- 
ening is always interrupted for T < Tg and the system 
gets stuck in one of those glassy states, from which it 
relaxes through a complex sequence of activated jumps. 
This explains the exponential increase of ri observed in 
Figl3] for T < 0.2. Once the system relaxes from the 
glassy state, it can be either directly equilibrate or decay 
first in a blocked state. 

In FigHH we show the typical behavior of relaxation 
function for g = 9 at a fixed temperature T < Tg and 
different system sizes. We see that the relaxation time 
(ri) is size independent for L > 200, which shows that 
the life time of the glassy states remains finite in the 
thermodynamic limit. 



D. Boundary conditions 

Finally, we analyzed the infiuence of the boundary con- 
ditions in the relaxation. To this end, we repeated some 
of the previous calculations using open boundary con- 
ditions. We found that the overall relaxation scenario 
found using periodic boundary conditions repeats quali- 
tatively for open ones. Moreover, the relaxation time as- 
sociated with striped configurations appears to be of the 
same order of magnitude of that corresponding to peri- 
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FIG. 14: (Color on-line) Relaxation function for ^ = 9, 
L = 200 and different temperatures around Tg. Tempera- 
tures increase from right to left. 
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FIG. 15: (Color on-line) Relaxation function for ^ = 9, 
T = 0.1 and different values of L. The inset shows a typical 
configuration of the glassy state associated with the plateau. 



odic boundary conditions. Although a more systematic 
study should be done to confirm that, it seems reason- 
able since the basic activated mechanisms here described 
should be still dominant in the case of open boundary 
conditions. In FigHHwe show an example of the equili- 
bration time probability distribution for q = 9 and some 
typical blocked spin configurations. In this case, Lif- 
shitz states are no longer composed only by hexagons 
for relatively small system sizes (due to the presence of 
the borders), but we see clearly the presence of stable 
three-colored vertex. Indeed, the observed configurations 
strongly resemble the blinking states reported in Ref^^^. 



IV. SUMMARY AND CONCLUSIONS 

The main conclusions of this work are summarized 
in the scheme of FigllTl After a quench from infinite 
temperature down to subcritical temperatures, the Potts 
model with single spin fiip kinetics and periodic bound- 
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FIG. 16: (Color on-line) Equilibration time probability dis- 
tribution P{r) for = 9, L = 200 and T = 0.2, using periodic 
(p.b.c) and open (o.b.c) boundary conditions. The inset im- 
ages show typical blocked spin configurations observed with 
o.p.c. 
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FIG. 17: Dynamical regimes in the long term relaxation of the 
q-state Potts model with ^ > 4, after a quench from infinite 
temperatures down to a sub-critical temperature T. 



ary conditions presents for g > 4 different relaxational 
regimes, determined by different crossover characteris- 
tic temperatures. Close enough to the critical tempera- 
ture, i.e., for Tn < T < Tc, relaxation is dominated by 
nucleation mechanisms. For intermediate temperatures 
T* <T <Tn the system crosses over into a simple coars- 
ening dominated regime where LAC law l(t) - t^/^ ^^^-yd^ 
until full equilibration, for most of the realizations of the 
stochastic noise. For lower temperatures Tg < T < T* 
the normal coarsening process is interrupted when the 
system gets stuck into highly symmetric blocked con- 
figurations, composed by macroscopic ferromagnetic do- 
mains, namely, striped and Lifshitz states. In those cases, 
the dynamic becomes activated with characteristic en- 
ergy barriers, which give rise to distinct time scales for 
the different process. 

Concerning the role of temperature in the relaxation 



through blocked states, we found that it has a double 
effect: at short time scales it enhances the probability of 
reaching them (which is zero at T = 0) and at long time 
scales it allows to escape from them through activation. 
At least for = 9, our simulations (for system sizes up 
to L = 500) suggest that the probability of reaching a 
blocked state at finite temperatures remains finite when 
L ^ oo. 

Striped states were previously found and characterized 
for the Ising model {q = 2) at very low temperatures. 
We found that their influence in the relaxation process is 
relevant for any value of q , even at relatively large values 
of T, but their occurrence probability is smaller for g > 3 
than in the Ising case. 

We found that the relaxation times associated with 
the blocked states present in general the finite size scal- 
ing behavior r ~ L^, where the exponent uo depends on 
T, taking values between 2 and 4. Such values of the 
exponent make the associated time scales several orders 
of magnitude larger than those associated with a normal 
coarsening process (which scale as L^) for large enough 
sizes, even at relatively large values of the temperature. 

Lifshitz prediction have been recently verified in the 
phase separation dynamics of diblock copolymers (Cahn- 
Hillard model), in a 2D hexagonal substrate^''. We ver- 
ified that Lifshitz prediction also holds for the q state 
Potts model with g > 3, even in a square lattice, if the 
system size is large enough. This strong finite size effects 
is probably due to the square symmetry of the lattice 
(large system sizes are required in order to the influence 
of the lattice to be faded out) and one should expect 
to be reduced in a lattice with three-fold symmetry (for 
instance, triangular). 

At very low temperatures T < Tg the system gets 
always trapped in glassy like metastable configurations 
whose life time is size-independent and diverge for T 
oo. After relaxation from the glassy state, the system 
can again gets trapped in a blocked state. Even when 
the glassy states do not dominate the relaxation at long 
enough time scales, a complete description of the relax- 
ation dynamics cannot exclude their existence and there- 
fore they deserve further investigations. 

Finally, we verified that the whole qualitative relax- 
ation scenario appears both for periodic and open bound- 
ary conditions, although the finite size scaling of the re- 
laxation times may differ in both cases. 
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